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Abstract 

Inspiraling compact binaries are promising sources of gravitational waves for ground and space- 
based laser interferometric detectors. The time-dependent signature of these sources in the detec- 
tors is a well-characterized function of a relatively small number of parameters; thus, the favored 
analysis technique makes use of matched filtering and maximum likelihood methods. As the param- 
eters that characterize the source model are varied so do the templates against which the detector 
data are compared in the matched filter. For small variations in the parameters, the output of 
the matched filter for the different templates are closely correlated. Current analysis methodol- 
ogy samples the matched filter output at parameter values chosen so that the correlation between 
successive samples is 97%. Correspondingly, with the additional information available with each 
successive template evaluation is, in a real sense, only 3% of that already provided by the nearby 
templates. The reason for such a dense coverage of parameter space is to minimize the chance that 
a real signal, near the detection threshold, will be missed by the parameter space sampling. Here 
we describe a straightforward and practical way of using interpolation to take advantage of the 
correlation between the matched filter output associated with nearby points in the parameter space 
to significantly reduce the number of matched filter evaluations without sacrificing the efficiency 
with which real signals are recognized. Because the computational cost of the analysis is driven 
almost exclusively the matched filter evaluations, a reduction in the number of templates evalua- 
tions translates directly into an increase in computational efficiency. Because the computational 
cost of the analysis is large, the increased efficiency translates also into an increase in the size of 
the parameter space that can be analyzed and, thus, the science that can be accomplished with the 
data. As a demonstration we compare the present "dense sampling" analysis methodology with 
our proposed "interpolation" methodology, restricted to one dimension of the multi-dimensional 
analysis problem. We find that the interpolated search reduces by 25% the number of filter evalua- 
tions required by the dense search with 97% correlation to achieve the same efficiency of detection 
for an expected false alarm probability. Generalized to the two dimensional space used in the 
computationally-limited current analyses this suggests a factor of two increase in computational 
efficiency; generalized to the full seven dimensional parameter space that characterizes the signal 
associated with an eccentric binary system of spinning neutron stars or black holes it suggests an 
order of magnitude increase in computational efficiency. 
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I. INTRODUCTION 



Inspiraling compact binaries of stellar mass neutron stars or black holes are among the 
most important gravitational wave sources accessible to the current generation of ground- 
based interferometric gravitational wave detectors 

mm 

. They are also very "clean" 

systems, in the sense that the gravitational wave signal arising from the inspiral depends 
only on general relativity (i.e., the structure of the binary components is unimportant) and 
can be calculated to great accuracy by the well-understood techniques of post-Newtonian 
perturbation theory 0,0, 3] • For these reasons, matched filtering and maximum likelihood 
techniques are well-suited for the detection and characterization of the signal from these 
systems jsl |^ and an implementation based on these methods is currently used in the 
analysis of data from the LIGO and GEO detectors [lo| . 

The gravitational wave signature of inspiral binary systems depends on a set of 15 pa- 
rameters that characterize the system (i.e., component masses, orbital energy and angular 
momentum at a given epoch, component spins, orientation relative to detector line of sight). 
To identify an incident signal using a matched filter requires the application of a fair sam- 
pling of filter "templates" , each defined by a unique choice of the parameters associated with 
the physical system. Current implementations of matched filtering used in the analysis of 
gravitational wave detector data involve a very dense sampling of the two-dimensional pa- 
rameter subspace corresponding to the binary component masses (intrinsic parameter space) 
and assuming zero eccentricity orbits and no body spins|3]: the templates are spaced so 
closely that the correlation between templates at neighboring points in the subspace is 97% 
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We refer to this as the "dense" search strategy. The rationale underlying the dense search 
strategy is to reduce the probability that a weak signal, characterized by parameters that fall 
between those sampled, will be missed by the sampling. Here we describe a straightforward 
and practical way of using interpolation to take advantage of the correlation between the 
matched filter output associated with nearby points in the parameter space to significantly 
reduce the number of matched filter evaluations without sacrificing the efficiency with which 
real signals are recognized. 

We are not the first to observe the significance of the high correlation between neighbor- 
ing templates nor to consider the opportunity for and advantages of interpolation as part 
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of the implementation of matched filtering for the analysis of binary inspiral signals. The 



significance of the high correlation as an indication that fewer templates shou 
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recover signals with the same efficiency, was first made in j3] • Croce et al 
the use of Cardinal interpolation with a truncated sine function to estimate the value of 
the matched filter output when the filter used corresponds to the actual parameters that 
describe the signal. They found a sampling of parameter space that would insure the inter- 
polated estimate would be no less than 97% of the maximum over a two dimensional intrinsic 
parameter space. Their sampling and interpolation reduced by a factor of 4, compared to 
the dense search, the number of templates required to search over a two dimensional intrin- 
sic parameter space. Here we find that we can achieve the same increase in efficiency per 
parameter space dimension with a simpler template spacing and a simpler and quicker to 
evaluate interpolation function. 

Other suggestions have been made for reducing the number of matched filter evalua- 
tions without sacrificing detection efficiency. One promising proposal involves a hierarchical 
search strategy, wherein a low-threshold trigger generated by the evaluation of the matched 
filters associated with a much coarser sampling of parameter space followed by (if necessary) 
a h i ghe r threshold evaluation matched filters over a much finer sampling of parameter space 
|. The interpolation proposal we make here is complementary in the sense 
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that it can be implemented together with the hierarchical strategies that have already been 
proposed to further improve the computational efficiency of binary inspiral analysis. While 
the gain in efficiency of the interpolated search over the dense search is approximately con- 
stant in the desired false alarm probability, the balance between the coarseness of the grids 
in the hierarchical steps, the number of hierarchical steps, and the gain in computational 
efficiency associated with the interpolation is not obvious and requires further study. 

The paper is organized as follows: In section |n] we describe the motivation behind our 
choice of interpolating function and the difference between our choice and the choice made 
in [l^ Q . In section IIIII we describe in detail the dense and interpolated search strate- 
gies, the two-dimensional template space used in current gravitational wave data analyses 
for inspiraling binary neutron stars, the one-dimensional restriction that we use here to 
compare the effectiveness of the interpolating search strategy, and (finally) compare the 
performance of the interpolated and dense search strategies by evaluating the sensitivity of 
each at fixed computational cost and the computational cost required by each to achieve the 
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same sensitivity. 



II. INTERPOLATING IN PARAMETER SPACE 

The Wiener matched filter W, corresponding to an expected signal characterized by r, 
is a scalar-valued function of the (vector-valued) instrument data d, noise power spectral 
density Sn- 

Wid\T)=W{T\Sn,d). (2.1) 

In our particular problem iy(d|T) is a continuous function of r and r corresponds to 
the parameters that characterize our binary system model: e.g., binary system component 
masses, orbital energy and angular momentum, component spins, etc. Given a data set d we 
wish to find an interpolating function W{t) and a set of points r^, in the space of possible 
signals such that 

Wk = W{Tk)=W{d\Sn,Tk). (2.2) 

There are, of course, an infinite number of continuous functions W{t) that take on the 
values Wk at the r^: the question is, how do we choose among them? 

Focus attention first on the case where r is a scalar x. One particular choice of interpolant 
W{d\Sn,x), which is especially important in the context of communication theory, is based 
on the Whittaker Cardinal function sine: 

oo 

C{x) = W.sinc^^, (2.3) 

A;=— oo 

where 

. , , sinvra; ,^ 
smc(x) = , (2.4) 

Tlx 

Xk = xo + kA. (2.5) 

Shannon showed that the Cardinal interpolation C{x) of W{d\Sn,x) is the unique 
interpolant W that (i) takes on the values Wk at the Xk, (u) has no singularities, and (Hi) 
and whose spectrum is limited to a bandwidth (2A)~^. Correspondingly, if W{d\Sn,x) is 
bandlimited in x and has the values Wk at the equidistant sampled points Xk then W{d\Sn,x) 
is equal to C{x). In the case where r is multi-dimensional the interpolation can be performed 
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separately on each index: e.g., in the case of two dimensions [i.e., r equal to (ri,r2)] 

oo 

C(t) = Wjksinc— (ti - tij) sine— (ts - rs,^) , (2.6) 

. , — ^1 ^2 

j,k=—oo 

where 

Tij = ri,o+jAi, (2.7) 
T2,k = r2,o + kA2 (2.8) 

and Ti 0, T2,o are constants. 

Cardinal inter pola tion using the Cardinal function sine forms the basis of the interpolation 
formula used in |l,4 ll^ . If W{d\Sn,T) is bandhmited and we choose our samples of W 
appropriately then we can do no better than using the Cardinal function to interpolate 
values of W between the samples. In our problem, however, iy(d|S'„, r) is not bandlimited 
and we do not have an infinite number of sample points Wk] correspondingly, the Cardinal 
function C(t) is at best an approximation to W{d\Sn,T). With that understanding the 
Cardinal interpolation C(t) is not preferred and we are led to seek other approximations to 

n 

W{d\Sn,T) that have favorable properties |30|. 

One possibility, chosen from approximation (as opposed to interpolation) theory, is the 
use of a Chebyshev polynomial expansion to approximate W{d\Sn,T). Without loss of 
generality consider a continuous function f{x) on [—1,1]. The Weierstrass Approximation 
Theorem states that for any e > we can find a polynomial P„ of order n such that 

max |/(x) -P„(x)| <e. (2.9) 

a;G[-l,l] 

The minimax polynomial approximation to W{d\Sn,x) is a natural candidate for the in- 
terpolation W{x). Unfortunately, finding the minimax polynomial is a very difficult process; 
nevertheless an excellent approximation to the minimax polynomial does exist. Define the 
error E{x\f,Pn) associated with the polynomial approximation P„(x) by 

E{x\f,Pn) = f{x)-Pn{x). (2.10) 

The Chebyshev Equioscillation Theorem states P* is the minimax polynomial if and 
only if there exist n + 2 points — 1 < a;o < a^i < ■ • ■ < Xn+i < 1 for which 

E{xk\f,P:) = {-l)'E, (2.11) 



where 



\E\ ^ max \E{x\f,Pn) 



(2.12) 



As a corollary, E{x\f, P*) vanishes for x G [—1, 1] at ri + 1 points x^, with Xk < x'^ < Xk+i- 
This result, together with the Mean Value Theorem, allows us to write the error term 
associated with the minimax polynomial P* as 



E{x\f,p: 

where ^ G [—1,1]. Correspondingly, 



(2.13) 



\E\ < max 

xe[-i,i] 



k=0 



k=0 



max 

Ce[-i,i] {n + l)\ 



(2.14) 



Focus attention on the order n + 1 polynomial 



Q 



n+l 



k=0 



X — Xi 



(2.15) 



This polynomial has leading coefficient unity. A unique property of the Chebyshev poly- 
nomial Tn+i is that, of all order n + l polynomials Qn+i with leading coefficient unity. 



max 

ze[-i,i] 



Tn{x) 



in— 1 



< max \Qn+i{x) 



(2.16) 



Additionally, T„+i(x) has exactly (n + 2) extrema on [—1,1], the value of |T„+i(x)| at these 
extrema is 1, and the extrema alternate in sign. Correspondingly, if the error term E{x\f, P*) 
associated with the minimax polynomial P* were polynomial — i.e., /^"'^'^■'(O were constant 
in equation 12. 131 so that i?(x|/, P*) was equal to Q*„ — then by the Equioscillation Theorem 
Qn^i would be equal to T^+i and the — where the error vanishes — would be the n+l 
roots of Tn+i- This suggests that we find the order n polynomial p* such that 



Pn(4) = VA; = 0...n 



(2.17) 



where, again, the x'f. are the roots of T^+i. The polynomial p* is a near minimax polynomial 
approximation to /(x). For this polynomial approximation Powell [2^ showed that, as long 
as /(x) is continuous on [—1, 1], 

■(^ + l/2)vr" 



-, . ^cheb . 
1 < < 

Co 



1 

1 H > tan 

n + 1 ^ 

fc=0 



2(n+ 1) 



(2.18) 



where 



eo = max |E(x|/,P„*)|, (2.19) 
echeb = max (2.20) 

a;G[-l,l] 



Powell also showed that grows slowly with n: in particular, 



2 

z/„ logn. (2.21) 

77 



Somewhat tighter bounds on z/„ can be placed when / is also differentiable 

As defined above, the near minimax polynomial is the interpolating polynomial that 
agrees with / at the n + 1 roots of T„+i. Alternatively, using several properties of Chebyshev 
polynomials, the Chebyshev interpolating polynomial can be expressed as a linear combina- 
tion of Chebyshev polynomials: 

n ^ 

K(x) = ^afcTfc(x)--ao, (2.22) 

fc=0 

where 



2 



a 



^ , (2.23) 



where, again, the are the n + 1 roots of T„+i. 



III. COMPARISON: DENSE AND INTERPOLATED SEARCH 

In this section we describe the dense and interpolating search strategy and compare their 
efficiency when applied to the problem of identifying the gravitational wave signature of 
coalescing neutron star systems in the LIGO detectors. 



A. Two Search Strategies 



The conventional search strategy used in the current analyses of LIGO, GEO and TAMA 



data (cf. 
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llll 12, [2^) begins with the placement of templates at discrete points Tk on 
the parameter space r. To choose the template locations we define the inner product of two 
signals g{t) and /i(t). 



{9,h) 



Snif) 



(3.1) 
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where g{f) is the Fourier transform of g and Sn is the detector noise power spectral density. 
Denoting by h{t\T) the signal characterized by r the match r(T_,-,Tfc) is 

T(t t)- {h{t\Tj),h{t,Tk)) .^2) 

" y/{h{t\T,),hit,T^)){hit\T,),hit,T,)y 

By construction |r| < 1. The templates locations are chosen so that consecutive templates 
in any of the directions tj have an overlap Tq, referred to as the "minimum match" and 
typically chosen to be 97%. 

With the templates placed, the dense search strategy proceeds: 

1. Evaluate the Wiener filter Vr(d|S'„, r^) at each of the template locations r^; 

2. Determine the template tj whose Wiener filter output is greatest; 

3. If the filter output at Tj exceeds the given threshold, report an event with the param- 
eters Tj. 

We refer to this as the dense search strategy. 

Following the discussion in section [H] we are in a position to describe an alternative 
strategy, which we refer to as the interpolated search strategy. First, fix the order of the 
interpolating polynomial. This determines the template locations on the parameter space 
T. Then 

1. Evaluate the Wiener filter W^(d|S'„, Tk) at each of the template locations r^; 

2. Form the interpolating polynomial from the iy(d|S'„, Tk)] 

3. Determine the location r' where the interpolating polynomial is maximized; 

4. Perform a final Wiener filter evaluation at r'; 

5. If the final evaluation exceeds the given threshold, report an event with the parameters 
t'. 

We illustrate the interpolated search strategy using Fig. [T| In Fig. [T] we use 37 interpo- 
lating search templates, that is, we sample the ambiguity function at 37 points in the tq 
space (the marked points on the dotted curve). We construct the interpolating function (the 
solid curve) and find its maximum by setting its derivative to zero. In order to avoid local 
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Demonstration of interpolated searcli using Chebysliev polynomials with 
37 templates over the range To=13-17sec and a test signal placed at Msec 
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FIG. 1: This figure demonstrates the interpolated search. The ambiguity function is sampled and 
reconstructed over the chosen parameter space of tq = 13 — 17 sec (only a part of the parameter 
space is shown) with the help of the Chebyshev interpolating polynomial. The approximate location 
of the peak of the interpolating function is first located and the zero of the derivative is obtained 
by applying successive approximations around the peak. Note that by placing a template at the 
maximum of the interpolating polynomial, the match has improved over the one obtained by simply 
evaluating the maximum of the interpolating polynomial. 



extrema, we first find the approximate location of the peak of the interpolating function and 
then find the zero of its derivative by successive approximation near the region of the peak. 
One can clearly see that the proper value of the ambiguity function at the maximum of the 
interpolating function is more than the maximum value of the interpolating function and 
this is what we gain by placing a template at the maximum of the interpolating function. 
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B. A one-dimensional parameter space for comparative studies 



We are interested in understanding the performance of the interpolated search strategy 
relative to the dense search strategy, which is currently used in the analysis of data from 
the LIGO, GEO and TAMA detectors 2- The current analyses focus on templates cor- 
responding to binaries with circular orbits and no component spins. The corresponding 
two-dimensional parameter space is spanned by the masses of the individual components. 
The templates vary most rapidly, however, along the axis spanned by the so-called chirp 
mass 

M := /i^/^M^/s, (3.3) 

where M is the system's total mass and fi its reduced mass. The linear density of templates 
needed by the dense search in the direction Om is approximately 100 times the linear density 
needed in the orthogonal direction. For the comparison we perform here we focus attention 
on the number of template evaluations needed for binaries with equal mass components 
that vary only in A4. We expect that the ratio of performance, measured as the number of 
templates required by the two search strategies to achieve the same search results, will be the 
same in the complementary dimension and in the other dimensions that will be introduced 
in future searches that accommodate component spins and orbital eccentricity. 



C. Templates 

The strain response of an interferometric gravitational wave detector to quadrupole for- 
mula approximation gravitational waves incident from an inspiraling binary neutron star 
system can be written 

h{t\ta,To) = ho [rtf{t -ta- To)Mf'\oS^{t - - Tq), (3.4a) 



where 



1/5 M ^^/^ 



-kM \ 256 TQ + ta — t 



$(t|t,,ro) := '^a + 2vr^ °rft/(t|4,ro) (3.4c) 

for t < ta + tq. Here ta is the moment when the instantaneous wave frequency is equal to 
fa and To is the elapsed time from that moment until (in this approximation) the system 
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coalesces, which is directly related to the system's chirp mass Ai: 

The elapsed time to coalescence tq is a useful surrogate for the chirp mass Ai: templates 
equispaced in tq have constant cross-correlation, independent of tq. Choosing fa equal to 
40 Hz, which is commonly taken as the lower-edge of the LIGO detector bandwidth at design 
sensitivity 2^, Tq ranges from approximately 43 s for a binary system consisting of two 1 Mq 
compact objects to 0.15 s for a binary consisting of two 30 M0 black holes. 

It is convenient to work with the Fourier transform of the strain response of the detector. 
For neutron star binaries in the LIGO or Virgo band the Fourier transform can be evaluated 
to an excellent approximation using the stationary phase approximation 

Hf) = J^r"' exp {i - 7r/4 + tq)]} , (3.6a) 

where 

*a = $(ta|ta,ro), (3.6b) 

vl>(/|4, ro) = 27r/4 + /„roy f -^J . (3.6c) 

The factor A/" is a constant amplitude. 

It is important to distinguish between the nature of the parameters that characterize the 
template. Changes in the parameter tq change the waveform shape: we term such param- 
eters dynamical parameters. On the other hand, parameters such as ta or $a translate the 
waveform, but do not alter its shape: we term these kinematical parameters. In our problem 
only the subspace of dynamical parameters needs to be spanned by discrete templates: the 
values of the kinematic parameters for the Wiener filter with the maximum output can be 
determined by other means. Correspondingly, at the level of approximation associated with 
the quadrupole formula the family of templates that must be evaluated is one dimensional. 



D. Dense search template placement 

There are many different ways of parameterizing the template space. Choosing tq as a 
dynamical variable has the advantage that r(ro, Tq) depends only on the difference tq — Tq; 
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Ambiguity function for Newtonian chirp using LIGO-i SRD with f|=40Hz, f^=512Hz 
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FIG. 2: Ambiguity function Ti plotted as a function of Atq. Horizontal lines are drawn for various 
matches; the dotted horizontal line is for a 97% match, which corresponds to an inter-template 
separation of ^ 30 ms. 

consequently, in the dense search templates are spaced uniformly in tq llE Ql- To 
determine that spacing we evaluate 



7^(Aro) =r(ro,ro + Aro), 



(3.7) 



where now F has been maximized over the kinematical parameters ta and (pa- This maxi- 
mization can be performed in a computationally efficient manner as shown in the literature 
We call Ti the dynamical ambiguity function or simply as the ambiguity function. 
It quantifies the fractional match between the template at tq and the signal at tq + Atq. 



Figure 
ments 



[shows Ti for power spectral density specified in the initial LIGO science require- 



J|. The requirement that 7i(Aro) is equal to a constant for any two consecutive 
templates determines the spacing Atq between templates that differ only in tq. For our 
example problem, which has just one dynamical parameter, the requirement that 7i(Aro) 
is 97% (the conventional choice) for neighboring templates leads to a template spacing Atq 
equal to 30 ms. 
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E. Interpolated search template placement 



In the dense search templates are equispaced in tq, with the spacing between adjacent 
templates — and thus the number of templates — chosen such that the dynamical ambiguity 
function takes on a specified value. When presented with data an event is signaled when 
the amplitude at one of these templates exceeds a threshold. 

In the interpolated search, on the other hand, the domain [r™™, t^^^] is mapped onto 
[—1,1] and the placement and number of templates is chosen to simplify the construction 
of the Chebyshev interpolating polynomial of the template output over this domain. When 
presented with data the maximum value of the Chebyshev interpolating polynomial is found 
and an event is signaled when the amplitude at that location exceeds a threshold. 

In the interpolated search our goal is to minimize the order of the interpolating polynomial 
(and, thus, the number of template evaluations) required for a given accuracy of interpo- 
lation. We have some control over this through the choice of mapping from [r™™, t^^^] to 
[—1,1]. The linear map 

r' = 2 ^° . - 1, (3.8) 

_max ^mm ' ^ ' 

is the most obvious such mapping. While have not made an exhaustive search of all pos- 
sible mappings; however, we have observed that better fits are possible with a lower-order 
polynomial when we use the mapping 



S = — cos 



^0 
TT 



(3.9) 



'0 '0 

Moreover, with this mapping, the roots of the Chebyshev polynomial are equi-spaced over 
the parameter range in tq. Once we have fixed the order n of the interpolating polynomial 
templates are placed at values of 6 that are roots of the T„+i(5) and the coefficients of the 
interpolating polynomial are found using equation 12.231 

In Fig. El we have plotted the match by placing normalized test signals (without noise) 
at regular intervals of tq. We see that the match is a (nearly) periodic function of tq, with 
the period equal to the template separation. This suggests that the detection probability is 
also periodic and this fact has been used in carrying out the simulations - the signals are 
injected within one such "period" in the parameter space. Moreover, one can see that with 
just 37+1 interpolated search templates one gets a minimal match of 0.97, whereas the 
dense search requires about 133 templates to achieve the same level of minimal match. This 
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for LIGO-I SRD with 37 templates 
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FIG. 3: Normalized test signals (without noise) were injected densely at regular intervals along the 
parameter space and the match obtained by the interpolated search strategy is plotted against tq. 
This figure demonstrates that the match is a (nearly) periodic function of tq with the period equal 
to the template separation. Moreover, with just 37 + 1 interpolated search templates the minimal 
match is 0.97. For the same minimal match 133 templates are needed for the usual dense search. 



amounts to a factor of 3.5 over the dense search and this is so for just one dimension. 



F. Comparison 

We are interested in two, related, comparisons: first, the relative "sensitivity" of a search 
carried-out with a fixed number of template evaluations using the dense search strategy 
and the interpolated search strategy; second, the number of template evaluations required 
by the interpolated search in order to achieve the same "sensitivity" as the dense search. 
To give meaning to the "sensitivity" of these two strategies we use the Receiver Operating 
Characteristic, or ROC. 
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The ROC is a plot of true positives as a function of the fraction of false positives for 
a binary classifier system as its discrimination threshold is varied. Both the dense search 
and the interpolated search are binary classifiers: i.e., they classify an interval of data d 
as including a signal or not including a signal. A true positive is a classification of d as 
including a signal when in fact it does; a false positive is a classification of d as including a 
signal when it does not. In both of the search strategies described here the discrimination 
threshold is matched filter output that must be exceeded for a data interval to be classified 
as including a signal. The false positive fraction is also known as the type II, or false alarm, 
error fraction and is denoted a. The fraction of true positives is also known as the detection 
efficiency e, which is one minus the type I, or false positive, error fraction (which is denoted 
/?). At fixed a a more sensitive search method has a greater e. The ROC associated with 
a search method no better than a toss of a (possibly loaded) coin is given by the diagonal 
a = e. 

Using numerical simulations we have evaluated a and e as a function of the detection 
threshold for both the interpolated search and the dense search, for different numbers of 
templates (dense search) and different interpolating polynomial order (interpolated search). 

To evaluate the false positive fraction a we generate a large number of data segments, 
each 2^^ samples long, and each consisting of Gaussian noise whose power spectrum (as- 
suming a 1024 Hz sample rate) is that specified as the initial LIGO science requirement 
j3| • (The Gaussian random numbers are themselves generated using the Mersenne Twister 
Pseudo Random Number Generator [3] and then filtered in the Fourier domain by scaling 
the Fourier components by the square root of the PSD.) For the purpose of this comparison 
we look for signals in the interval tq G [13 s, 17 s]. Both the dense and interpolated search 
methods are applied to this data. The ratio of the number of events signaled to the number 
of data segments examined as a function of the threshold is a for that threshold. Ap- 
proximately 50,000 realizations of detector noise are used to evaluate a, which gives reliable 
results for a greater than approximately 10~^. 

To compute e, the true positive fraction, we proceed in a similar fashion. Now, however, 
with each noise instantiation we add a signal, with tq drawn uniformly and randomly from 
the interval covered by the search: i.e., tq G [13 s, 17s]. In almost all cases 50, 000 realizations 
of detector noise plus signal are used to evaluate the efficiency, which gives rehable results for 
efficiencies greater than approximately 10^^. However, for the fiat search with 40 templates 
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and the interpolated search with 30 templates,we have used 400, 000 reahzations. The larger 
number of realizations in these cases results in smoother curves. 

The top panel of Fig.|3]shows the variation of a for both methods using 40 templates: i.e., 
a 100 ms template spacing for the dense search and an order 39 interpolating polynomial 
in 5 (cf. equation 13. 9|) . For any threshold a is always greater for the dense search than 
for the interpolated search; similarly, as shown in the center panel of Fig. |3J for any given 
threshold the efficiency e is always greater for the interpolated search than for the dense 
search. Finally, the bottom panel of Fig. 0] shows the ROC for a 40 template dense search 
and an order 39 interpolated search, both of which involve 40 template evaluations to decide 
if a signal has been detected. Comparing both ROCs it is clear that the interpolated search 
is more sensitive at any given a then the dense search. This is always true: i.e., for a fixed 
number of template evaluations the interpolated search will always have a better efficiency 
at a given a than the dense search, though as the number of templates grows large the 
fractional difference in sensitivity will decrease. 

Figure El and table I addresses the second of our two questions: the number of templates 
evaluations required of an interpolated search to have the same sensitivity as a dense search. 
Figure El shows the ROCs for dense searches using 140 and 160 templates, together with 
the ROCs for interpolated searches using 120 and 100 templates. The interpolated search 
with and order 120 interpolating polynomial is clearly as sensitive as a dense search with 
160 templates, and an interpolating search with an order 100 polynomial is as sensitive as a 
dense search with 140 templates. Table I shows similar pairings of the number of templates 
in a dense search and the number of templates in an interpolating search necessary to achieve 
the same sensitivity. 

IV. CONCLUSION 

We have shown that the use of near-minimax interpolating polynomials to fit the output 
of matched filters to the filter parameter values can greatly improved the sensitivity of a 
matched-filter based search for gravitational waves from compact binary inspiral. Using 
such a polynomial to find the parameters of the signal template leading to the best match 
we can reduce the computational cost of a search over a two dimensional parameter space 
by a factor of two compared to the methods currently in use, without any loss of sensitivity 
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Comparison between false alarm probability vs. threshold curves 
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Comparison between false dismissal probability vs. threshold curves 
for injected SNR=8 
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FIG. 4: The variation of the false and true positive fractions, a and e with threshold r] for the 
dense and interpolated search methods, each making use of 40 template evaluations. The top panel 
shows the false positive fraction. Note how the false positive falls much sooner for the interpolated 
search than for the dense search. The bottom panel shows e when a signal of amplitude signal to 
noise 8 is present in the range tq S [13 s, 17 s]. Note how the e is always greater for the interpolated 
search than for the dense search. For the same computational cost (determined by the number of 
template evaluations) the interpolated search will always perform better than the dense search. 
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TABLE I: Number of template evaluations required to obtain the same efficiency at a false alarm 
fraction of 10~^ in a dense search and an interpolated search. Note how the interpolated search is 
computational more efficient for the same sensitivity. 



# templates 


e at a = 10"^ 


Dense 


Inter p. 


40 


31 


0.859 


50 


41 


0.890 


60 


49 


0.905 


80 


64 


0.919 


100 


89 


0.924 


140 


105 


0.927 


160 


115 


0.929 



or discriminating power. This factor of two becomes a factor of ten when the search is over 
the seven dimensional parameter space that includes not only the masses but also the spins 
of the binary components 



28| . This savings in computational cost is estimated under the 
assumption, which we believe well-founded, that we will obtain the same savings when the 
interpolation is extended to additional dimensions. 

The use of near-minimax interpolation should be considered as part of a larger strategy 
that employs a multi-grid approach to determine whether a signal has been observed and, if 
so, the parameters that characterize it. Since the major contribution to the computational 
cost of a multi-grid search is thought to arise in the initial stage of the search the gain in 
computational efficiency — and, correspondingly, the size of the parameter space that can 
studied with fixed computational resources — could be substantial. 
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FA vs FD probability curves for SNR=8 with coarse bank of templates 
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FA vs FD probability curves for SNR=8 with fine bank of templates 
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FIG. 5: ROC curves for dense searches (solid curves) and interpolating searches (dashed curves). 
For a given number of templates, the solid curves are 'lower' - less false dismissal probability for 
the same false alarm - than the dotted curves in the regime of low false alarm showing that the 
interpolated search performs better than the dense search for low false alarm. The bottom panel 
shows analogous plot for high minimal match (fine bank) ~ .98. Here the performance of the dense 
search with 160 and 140 templates is comparable to that of the interpolated search with 120 and 
100 templates respectively. 
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